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Abstract: We derive the fundamental equations describing nonlinear 
propagation in multi-mode fibers in the presence of random mode coupling 
within quasi-degenerate groups of modes. Our result generalizes the 
Manakov equation describing mode coupling between polarizations in 
single-mode fibers. Nonlinear compensation of the modal dispersion is 
predicted and tested via computer simulations. 
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1. Introduction 



Optical spatial division multiplexing in multi-mode optical fibers (MMFs) has been recently 
attracting enormous attention as it is considered to be a promising solution for the current ca- 
pacity crunch. Optical nonlinearity is one of the most fundamental aspects of fiber propagation 
and its effect on the ability to transmit information is of utmost importance. In a recent paper 
we have considered the case of nonlinear propagation in multi-mode fibers where all excited 
modes experience strong random mode coupling [T|. In realistic fiber-communications scenar- 
ios, it often happens that while some modes have similar wavenumbers and therefore experience 
strong random coupling, their coupling to other modes (whose wavenumbers are significantly 
different) is much weaker. The case of the weakly guiding step-index optical fiber [12 is a typ- 
ical example for this situation as the linearly polarized LP m „ modes contain 2-fold and 4-fold 
degeneracies. In this paper we show that in this framework nonlinear multi-mode propagation 
is accurately described in terms of coupled generalized Manakov equations (ME) Q and give 
the expression for the coupling coefficients. The advantage of this description over that based 
on standard coupled nonlinear Schrodinger equations (NLSE) is that it significantly reduces 
the number of relevant nonlinear parameters and hence simplifies the entire description of non- 
linear propagation. It also forms an excellent starting point for analytical studies of nonlinear 
effects in fiber-optic transmission |4|. In order to validate the accuracy of the coupled ME, we 
show that the solitary solutions that are predicted on their basis can indeed be observed in direct 
simulations of the complete model which is based on the coupled NLSE. 

2. Analysis and results 

We start from the most general form of the coupled NLSE for MMFs |5J and analytically intro- 
duce the effect of random mode coupling within groups of degenerate modes. The assumption 
under which our theory works is that the correlation length characterizing the effect of random 
mode coupling is shorter [ 1 ] than the length-scale of the nonlinear interaction between pulses. 
This regime is strictly satisfied in the case of the single-mode fiber supporting two polarization 
modes (6) and since the mechanism of mode coupling is caused by perturbations which are of 
the same nature, extension to larger mode numbers (in the absence of explicit experimental data 
on the nature of perturbations) is plausible. This assumption is also consistent with reported ex- 
perimental results 00. Our formalism allows a seamless generalization of all results known 
in the context of nonlinear propagation in single-mode fibers to the multi-mode case. 

We consider a fiber with N spatial modes, such that the total number of scalar modes in- 
cluding polarizations is 2N. Defining a 2A^-dimensional electric field vector E such that its 
components are the complex amplitudes of the various modes, the coupled NLSE describing 
field propagation through the MMF is given by 

d -f = iB(°)E~BW^-i^^+i 7 £C Plkm E; 2 E k E m ej, (1) 

jhkm 

where the dimensionless constants Cjhkm depend on the details of the spatial mode profiles. 
We find it convenient to introduce y = (iXsnij cA e g, which is identical to the usual nonlinearity 
coefficient appearing in the scalar NLSE of a single mode fiber, where n% is the Kerr coefficient 
of glass, c is the speed of light in vacuum, and A e g is the effective area of the fundamental 
mode at central frequency cOq. Here, BW(z), i = 0,1,2 are 2N x 2N Hermitian matrices of 

components (9) (z) and we use ej with j = 1 , . . . , 2N to denote the set of complex orthogonal 

unit vectors used to represent the electric field, namely E = J^jEjij. We now assume that the 
modes are divided between two groups denoted by a and b. The two groups are characterized 
by distinctly different propagation constants, but the modes within each of the groups are nearly 



degenerate. We will denote the number of degenerate modes in groups a and b by N a and Nb, 
respectively, such that N a + N b = N. We assume the existence of only two such groups only 
for the simplicity of presentation, since the final results can be readily expanded such that any 
number M of groups of degenerate modes can be accounted for. We denote by j3 ( , and by f5 b 
the wavenumbers of the two groups of modes, and assume that the difference between them 
is large enough to ensure that mode coupling between the groups is negligible on the scale of 
the nonlinear length. In the case of the LPoi and LPn mode groups of a step index fiber, this 
property has been verified both computationally |9l and experimentally Q HD . The coupled 
generalized multi-component [10| ME, which are the main result of this Letter, are given by 
ED: 

^ = i^-jS^-i^^ + iy^l^ + ^l 2 )^ (2) 
= iP b E h -l3^-i^^+ir[K ba \E a \ 2 + K bb \E b \ 2 )E b , (3) 

where E a and E b are the vectors describing the modes in groups a and b, respectively, in a 
reference frame that accommodates for the unitary evolution induced by the linear coupling 
within each group, and their respective number of components is 2N a and 2N b . The quantities 
b and f}" b are the group velocity and the group velocity dispersion terms of the two groups 
of modes, and the parameters K uv are generalized self-phase modulation and cross-phase mod- 
ulation (XPM) coefficients given by 

where 5,-,, is the Kronecker delta function. In @, each of the indices u and v takes the values a 
and b depending on which of the coefficients in Eqs. (f2]i and © is evaluated. For example, if 
u = a and v = b, then in the summation the index j runs over all values corresponding to the 
modes in group a and the index h runs over all the indices corresponding to the modes in group 
b. The indices k and m in the summation span all the modes. Because of the symmetry of the 
coefficients (5] C jhkm = C hjkm we have K uv = K vu = /X. 

The coupled ME (ffl and ® have a significant complexity advantage over the generic NLSE 
Q), which is of relevance both in the numerical and analytical contexts. In the numerical con- 
text, we note that in a fiber supporting M groups of degenerate modes, there are M self -phase 
modulation coefficients Kjj and M(M — l)/2 cross-phase modulation coefficients, one for every 
pair of groups. The overall number of nonlinear coefficients thus scales quadratically with the 
number of groups. This is to be compared with the number of coefficients Cjhkm appearing in 
Eq. (Q]i, which scales as the fourth power of the total mode number N. In the analytical con- 
text, we note that while Eq. (fTJ in its generality does not display any simple qualitative feature 
(in general, solitary pulses are not supported), we will show in what follows that Eqs. (f2]i and 
(0 have the characteristic behavior of two coupled nonlinear equations that support soliton 
solutions. 

The derivation of Eqs. (ffl and (0 starts from the notion that the nonlinear interaction between 
the two groups of modes can be averaged with respect to the random linear evolution of the 
fields within each group. That is because in the regime of strong coupling the linear evolution 
within the groups is significantly faster than the nonlinear evolution. As the statistics of the 
linear evolution are isotropic, the averaged nonlinear coupling is bound to be isotropic as well, 
and the nonlinear coupling between the groups can only depend on the fields amplitudes, not 
on their orientations. In this case, the nonlinear term in the averaged equations representing 



the evolution of the field in mode u must have the form iyyK uv \E v \ 2 + K uu \E u \ 2 j E u , where 
the first and second terms correspond to the effects of cross-phase and self -phase modulation, 
respectively. Consistency with the last term of Eq. ([T]i implies the validity of the following 
substitution 

52 CjhkmElE k E m ej i — > Kiv\E v | E u . (5) 

k,m j€uhEv 

Scalar multiplication of the right-hand-side of (0 by E u is equivalent to scalar multiplication 
of the left-hand-side by the extension of E u to the 2N dimensional vector representation used 
in Eq. (fTJ. Averaging with respect to the orientations of E u and E v yields 

p*p,p p* 

Kiv = 2-2-2- Cjhkm , 2 .0 |2 ; (6) 

k,mjeuh£v \^u\ \£-v\ 

where the symbol £ denotes statistical averaging. The derivation of Eq. from Eq. © is 
presented in the appendix. Equations ©-(IUi have been derived in |12| for the special case 
N a = Nfr = 1 using a different approach. This special case represents a situation in which none 
of the spatial modes couple with each other, so that strong coupling can only be attributed to 
polarizations. 

3. Numerical validation 

In order to validate Eqs.(f2]i-(|4]i, one would have to demonstrate that for any arbitrary excita- 
tion, propagation according to these equations leads to the same field evolution as propaga- 
tion according to Eq. (Q~|). In principle this can be done with the kind of excitation commonly 
encountered in telecom applications, but a far more convenient and elegant choice would be 
to take advantage of the fact that Eqs. ©-(HJi can be shown to have analytical solitary solu- 
tions with predictable properties. In the case of solitary pulses, nonlinearity plays a much more 
pronounced role than it does in most typical communications scenarios. Therefore, by demon- 
strating that the same pulses are observed in the numerical solution of Eq. (Q]), the accuracy 
of Eqs. ©-(HI can be verified. In the numerical studies discussed in what follows we assume 
a weakly guiding step-index fiber supporting the LPoi and LPn mode groups. Since mode 
LP] i is two-fold degenerate, this scenario is represented in our notation by setting N a = 1 and 
Nb = 2, where a and b refer to modes LPoi and LPn, respectively. When only one of the groups 
is excited (in a particular mode and polarization), it can be shown by direct substitution that a 
hyperbolic-secant pulse of the form 

E u = a„sech ^ tT^ ) eXp (~ ia)ut } ( - 7 - ) 

is an exact solitary solution of Eqs. where a u , T u , x u and co u are parameters representing 

the amplitude, the central position, the temporal width, and the center frequency — defined as the 
offset from the carrier frequency — of the optical pulse and where a u and X u satisfy the relation 
T u fl„ = -y/ 1/3// 1 /(/£„„) and co u = 0, with u = a or u = b. In the case where both mode-groups 
are excited, it can be shown that solitary solutions still exist and can be very well approximated 
by pulses whose functional form is identical to @, with the parameters a u , T u , x u and co u (for 
u = a,b) being analytically expressible |4l. 

In our simulations, we excited one polarization mode in each of the LPoi an d LPii mode 
groups with a hyperbolic secant pulse as in (0. The parameters a u , T u and x u were set to their 
predicted steady state values, but the center frequencies at the launch position were both set to 
zero (while their steady state values would be different from each other H). The pulses were 




Fig. 1. (a) Three-dimensional plot of the power envelope of the two modes obtained from 
Eq. Q}. The fact that solitary evolution is demonstrated validates the coupled ME, on the 
basis of which solitary evolution was predicted, (b) Timing difference T a — T/,, center fre- 
quencies (denoted by ft),, and a>b), and pulse-widths t„ and % as a function of the propaga- 
tion distance z. The black-dashed curves represent the solution of Eq. Q} and the solid-red 
by the solution of the coupled Manakov equations (O-©. These two curves are hardly 
differentiable in the resolution of the figure. The horizontal dashed-dotted lines show the 
analytical steady state values obtained in |4). 

propagated based on the full NLSE model (Q]) with the expectation that they should converge 
to the solitary solution. We assumed a core radius of 7.5/im, a core refractive index of 1.46, a 
refractive index step of 9.7 x 10~ 3 (similar to lfl3l ). dispersion coefficients ft" = — 25ps 2 /km 
and fi' b ' = — 29ps 2 /km, and a differential group delay between the two groups of j3,' — fi' b = 
0.2ps/km. The nonlinear coefficient in our computations was «2 = 2.6 x 10 _20 m 2 W _1 and 
the effective area for the fundamental mode was A e ff = 126/im 2 , corresponding to y ~ 0.835 
W _1 km _1 , K aa ~ 0.89, K bb ~ 0.76, ~ 0.895. A fundamental difficulty in the simulations 
is to accommodate the very large wavenumber difference between the two groups, which is 
estimated to be j5 a — Pb — 9.153 x 10 3 m _1 . This value corresponds to a beat-length smaller 
than a millimeter and since the overall length of the simulated fiber should be in the tens or 
hundreds of kilometers range, simulations become prohibitively inefficient. In order to bypass 
this difficulty, the wavenumber that we used in the simulations was scaled down by a factor 
of 10,000, resulting in a beat-length of the order of 10m. The integration-step in the split-step 
solution of the coupled NLSE was correspondingly limited not to exceed 0. lm. 

Figure [TJ a) shows the power envelope of the two pulses propagating in the LPoi and LPn 
modes (denotes as \E a \ 2 and |£i,| 2 , respectively). These plots were obtained by solving Eq. 
(Q]l numerically. The solitary evolution of the two pulses, which is predicted on the basis of 
the coupled ME ©-(O ID, validates the accuracy of the latter. In Fig. [TJb) we show several 
attributes of the pulses propagating in the two groups of modes. The initial temporal widths of 
the two pulses used in the simulation were x a = 9.3ps and z b = lO.lps, and the pulse energies 
in the LPoi an d LPn modes were 5pJ and 2.5pJ, respectively. In each of the plots in Fig. 
[TJb) the dashed and the solid curves represent the results of the coupled NLSE ([JJ) and of the 
coupled ME ©-(O respectively. As is evident from the figure, the difference between those 



two curves is hardly noticeable indicating the accuracy of the coupled Manakov model. The 
fiat dotted lines represent the parameters of the steady state solution described above to which 
the numerical solutions are clearly seen to converge. In the steady state solution, the initial 
difference of group velocity of the two groups of modes is compensated by a self-generated 
frequency shift of the two solitons [4|. This dynamics is identical to soliton trapping observed 
in single-mode birefringent fibers 1 14 1. 



4. Conclusions 

To conclude, we have derived a set of coupled ME describing the nonlinear evolution of the 
electric field in a MMF in the presence of random mode coupling. These equations highlight 
that the nonlinear evolution in MMFs in the presence of random mode coupling is characterized 
by a small number of parameters that scales with the square of the number of non-degenerate 
groups of modes. This is in striking contrast with the the number of parameters required to de- 
scribe the nonlinear interaction in the absence of random mode coupling, which scales with the 
total number of modes to the fourth power. In order to demonstrate the validity of the coupled 
ME, we showed that an explicit numerical solution of the full coupled NLSE (Q~|i produces the 
solitary solutions predicted by the coupled Manakov model. 

Appendix: Equivalence of Eqs. © and ©. 

Our goal is to perform the average appearing in Eq. © with respect to the orientations of the 
filed vectors E u and E v , which can be assumed to have uniform distributions. To this end we 
introduce auxiliary random vectors X u and X v , with 2N U and 2N V complex, statistically inde- 
pendent components. The real and imaginary parts of each component Xj of X u and X v are sta- 
tistically independent standard Gaussian variables having zero-mean and unit variance. Using 
known properties of Gaussian vectors, it can be shown that the term multiplying the coefficients 



Cjhkm in Eq. © is equal to Q hkmj = S 



I ! \Xv 



/(\Xu\ 2 \X v \ 2 ) 



Multiplying both 



Y* Y Y Y* 

h * m < j 

sides by |X„| 2 |X V ,| 2 and performing another average (with respect to the square modulus \X U \ 2 
and |X V | 2 ), we find that 

S 

■ 

(8) 
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\t\ 2 \x^ 



(2N U )(2N V + 8 VU ) : 



where we have made use of standard properties of Gaussian variables. 



